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Abstract. In a system of coupled harmonic oscillators, the interaction can be represented 
by a real, symmetric and positive definite interaction matrix. The quantization of a Hamil- 
tonian describing such a system has been done in the canonical case. In this paper, we take 
a more general approach and look at the system as a Wigner quantum system. Hereby, 
one does not assume the canonical commutation relations, but instead one just requires the 
compatibility between the Hamilton and Hcisenbcrg equations. Solutions of this problem 
are related to the Lie superalgebras fl^lln) and osp(l|2n). We determine the spectrum of 
the considered Hamiltonian in specific representations of these Lie superalgebras and discuss 
the results in detail. Wc also make the connection with the well-known canonical case. 
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1 Introduction 

In quantum as well as in classical mechanics, the harmonic oscillator is one of the most popular 
examples to describe harmonic movement of a particle. Their numerous applications and their 
analytical solvability as quantum systems explain why harmonic oscillator models are thoroughly 
investigated. Systems of interacting harmonic oscillators are among these well-known models. 
A system of n one-dimensional harmonic oscillators interacting with each other can be described 
in its most general form as [T| 

H = r'^Vr. 



In this equation, r ' is the vector {p\ 



ql), with qj. and Pr respectively the position 



and momentum operator of the oscillator at location r, and y is a positive definite matrix 
describing the coupling in position and momentum coordinates. In the present paper, we will 
assume that there is no coupling involving the momentum operators. Following this approach, 
we can write the Hamiltonian in the following manner: 



\Pr. 



(^1 



The matrix A is called the interaction matrix, and it is assumed to be real, symmetric and 
positive definite. In order to connect the physical context of harmonic oscillators coupled by 
springs obeying Hooke's law to this Hamiltonian, we can rewrite A as uj'^I + cM. All oscillators 
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then have mass m and natural frequency lo, and the couphng constant is cahed c (c > 0). The 
nx n identity matrix is denoted / and M is a general real and symmetric matrix. This notation 
is solely introduced to be able to interpret the system physically. The essential mathematics can 
and will be done using the more general notation A. 

Much interest lies in the spectrum of the Hamiltonian ([T|), since this yields all the possible 
values that might arise when measuring the energy of the system. In the standard approach for 
determining this spectrum, one imposes the canonical commutation relations (OCRs): 



This has been done for several types of interaction matrices in [2]. However, there is a more 
general approach to tackle this problem. Eugene Wigner was the first to realize that one does 
not need to assume the OCRs in order to find operators that satisfy Hamilton's equations (in 
operator form) and the equations of Heisenberg simultaneously. Instead, imposing that these 
equations are equivalent as operator equations results in a set of compatibility conditions (CCs). 
In a standard quantum system, these CCs are naturally satisfied as a consequence of the CCRs. 
Wigner on the other hand questioned the fact that the relations ([2]) can be derived from the 
compatibility conditions [3]. His discovery that this was not the case for a single harmonic 
oscillator, resulted in the first Wigner quantum system [3] . 

In the present paper, we consider a quantum system of n coupled one-dimensional harmonic 
oscillators, and we treat this as a Wigner quantum system. We will demonstrate how the analysis 
of the spectrum of this system is connected to the Lie superalgebras osp(l|2n) and 0((l|n). The 
Wigner quantization of this string of harmonic oscillators deviates from its standard quantum 
mechanical counterpart by a certain parameter, corresponding to the parameter characterizing 
unitary irreducible representations of these Lie superalgebras. The representations that we will 
consider in this paper are both characterized by a parameter p. For 05p(l|2?i) one finds back 
the canonical case by choosing p = 1. 

It is interesting to note that Wigner in his original paper [3] already found a connection 
with osp(l|2), although he was unaware of this fact. His compatibility conditions were exactly 
the defining triple relations of this Lie superalgebra, or equivalently, the triple relations of the 
paraboson algebra. Green [5] later generalized these results and was the first to write down the 
paraboson relations, connected to 05p(l|2n), explicitly. 

Other scientists have been inspired by Wigner's approach, which has resulted in the research 
of several different Wigner quantum systems O [71 [8]. Amongst them, a system of coupled 
harmonic oscillators with periodic boundary conditions has been studied in [9], where solutions 
for the position and momentum operators are found in terms of generators of the Lie super- 
algebra 0[(l|n). Analysis of the properties of this quantum system has been done in a Fock 
type representation space of g[(l|n) and the authors found a discrete and finite spectrum of 
the coordinate and energy operators. Another quantum system consisting of coupled harmonic 
oscillaters with a fixed wall boundary condition has been the subject of investigation in |10j . 
Here the authors present solutions in another class of representation spaces of 0l(l|n), called the 
ladder representations. These systems of coupled harmonic oscillators, however, correspond to 
specific interaction matrices A. We wish to extend the performed analysis to a general inter- 
action matrix and investigate properties of the system in specific representations of g[(l|n) and 
osp(l|2n). 

This paper contains two specific examples of interaction matrices that will get much attention 
throughout the manuscript. The first system is a system with constant interaction, and its 
Hamiltonian is 



[qr,qs] = 0, 



[Pr,Ps] = 0, 



[qr,Ps] = ihSrs- 



(2) 
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with % = Qn+i = (fixed- wall boundary conditions). Much is known already about this system. 
We refer to [TO] for more information. The second example is the Hamiltonian with so-called 
Krawtchouk interaction. This Hamiltonian depends on a parameter p, and for p = \ one ends 
up with its simplest form: 

1 .0 / 2 c(n — 1) \ .0 cm v-^ /—. r. ^ ^ 

^^2^ Z^Pr + ^Z^(^ 2 J ~ ~ ^ ^^('^ - r){qrqr+i + qr+iqr)- 

r=l r=l ^ ^ r=l 

Notice that we cannot say that qr-qr+i + qr+iqr- = 2qr'?r+i because we make no assumptions 
about the commutation relations between the position and momentum operators. In Section [3] 
and onwards we give the interaction matrices that belong to the considered systems and explain 
why these examples are interesting. 

In Section [2] we translate the problem into a different form, which is connected to Lie super- 
algebras. This analysis works for an interaction matrix in its most general form. We supply the 
reader with concrete examples of interaction matrices in Section [3l The given examples handle 
linear chains of coupled oscillators and analytically solvable interaction matrices. This means 
that the procedure from the previous section works analytically, not only in a numerical way. 
Section [J] restates the problem in terms of Lie superalgebra generators. Hereafter, one can try 
to determine the actual spectrum in specific representations of the Lie superalgebras g[(l|n) 
and osp(l|2n). This is done in Section [5l where the reader is supplied with general formulae to 
determine the spectrum and a detailed analysis including some plots. Finally, we go back to the 
known canonical quantization and find connections with the previous results. 



2 The Wigner quantization procedure 

The Hamiltonian ([T]) can also be written as 

1 " ?Ti " 

^ = ^Y.P^r+'^^arsqrqs. (3) 

r=l r,s=l 

where we have introduced the notation a^s for the element on position (r, s) of the interaction 
matrix A. We also assume that that the position and momentum operators are self-adjoint, that 
is ql = (jr and pi = Pr- Instead of imposing the canonical commutation relations ([2]), we just 
require the equivalence of Hamilton's equations (in operator form) 

, _ dH , _dH 

Pr — 77^ J qr — ^ . j 

Oqr OPr 

and the Heisenberg equations 

Pr = j\H,Pr], qr = -[H,qr]. 

The resulting compatibility conditions, applied for the Hamiltonian ([3]) become 

ih " 

[H,qr] = Pr, [H,Pr] = ihmS^ Qrsqs, (4) 

m ^-^ 

with r = 1,2, ... ,n. We are now looking for operator solutions for (jr and pr satisfying the 
compatibility conditions (UD, with Hamiltonian Q. Since the interaction matrix is real and 
symmetric, we can apply the spectral theorem to A and write 

A = UDU^. 
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In this identity, D is a diagonal matrix with the real and positive eigenvalues fxj (j = 1, . . . ,n) 
of A as diagonal elements. The columns of the matrix U are the real and orthonormal eigenvec- 
tors of A, is the transpose of U. Of course, U is an orthonormal matrix. In other words it 
satisfies 



I 



U^U. 



The matrix U can be used to introduce the normal coordinates and momenta. These new 
operators are defined by 



u 



T 



(Pi 



\Pn 



(5) 



The operators Pj and Qj are self-adjoint and they do not satisfy the OCRs, just like the opera- 
tors Pr and Qr- In function of the normal coordinates and momenta, the Hamiltonian ([T]) can be 
rewritten as 



H 



1 

2m 



[PI 



/A 



Qk)D 




or, more explicitly in summation form 

^ n n 



(6) 



This Hamiltonian shows that only the eigenvalues of the interaction matrix A are of essence 
when the system is rewritten in function of the new operators Pj and Qj . Looking for operator 
solutions of the Hamiltonian ([6]), one should not forget to take the compatibility conditions (jU 
into account. For the normal coordinates and momenta, these translate into 



ih 



[H, Qj] = Pj, [H, Pj] = ihmfijQj 



m 



(7) 



This can be obtained by substituting the transformations ([5]) in the compatibility conditions ^ . 

It turns out that we will be able to find solutions for Qj and Pj satisfying the CC's ([7]) and 
the Hamiltonian in equation ([6]) in terms of Lie superalgebra generators. The easiest way to 
establish such a result, is to introduce linear combinations of the unknown operators Qj and Pj 
as follows: 



qij 



2h 



^J^hnuJJTj 



--P, 



(8) 



In terms of the operators , which satisfy the adjointness relations (a^- ) 
nian ([6]) can be rewritten as 



aj, the Hamilto- 



n t I 



n ^ I — - 



(9) 



Again, we need to have the compatibility conditions in terms of the newly introduced operators. 
These follow from d?]) and are 



1,2, ... ,n. 



(10) 



Thus we have: 
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Theorem 1. The Wigner quantization of the system ^ has been reduced to the problem of 
finding 2n operators (j = 1, . . . acting in a certain Hilbert space. These operators must 
satisfy (o^)^ = a J and 

n 

^ [^//I7{a+,oj},a^] = ±2^//i^a^, fc = l,2, ...,n. (11) 

The Wigner quantization procedure is reversible, so that the knowledge of the operators allows 
us to reconstruct the observables Pr and The Hamiltonian is given by equation ([9]). 

Equation (jlip is equivalent to a quantum system describing an n-dimensional non-isotropic 
oscillator [11^ Section 2]. For such systems, it is known that solutions in terms of Lie superalgebra 
generators exist [11]. Some specific solutions are related to the Lie superalgebras 05p(l|2?i) 
and g[(l|n), but not all solutions are known for n > 1. We will focus on these two solutions and 
investigate the spectrum of our system in representations of these Lie superalgebras. However, 
before moving on to this analysis, we will give some explicit examples of interaction matrices. 



3 Some interaction matrices 



The Wigner quantization procedure only requires the spectral decomposition of the interaction 
matrix. Since we assume that A is real and symmetric, the spectral theorem is always appli- 
cable. Hence, for a given interaction matrix, the Wigner quantization procedure always works 
as above. However, the explicit spectral decomposition has to be calculated numerically. For 
some matrices, we have analytically closed expressions for the eigenvalues and eigenvectors of A. 
Hamiltonians described by such interaction matrices are called analytically solvable, and they 
were thoroughly investigated in [2]. We will consider two examples of such analytically solvable 
interaction matrices. The interaction of the resulting systems will be referred to as constant 
interaction and Krawtchouk interaction. 

Constant interaction jl2|, I13^ I14j . The interaction matrix ^cst = ^"^I + cMcst, with Mcst 
the n X n matrix given by 



est 



/ 2 



V 



-1 



-1 



-1 



\ 



is an analytically solvable interaction matrix. It represents a linear chain of coupled identical 
oscillators with constant interaction throughout the chain. The corresponding Hamiltonian is 



He 



r=l 



Pr 

2m 



+ 



muj 



Ecm , ^ 

r=0 



2r+l 



with qo = qn+i = 0. The interaction between the oscillators is of a nearest-neighbour type, 
which is a direct consequence of the tridiagonal form of the interaction matrix. The eigenvalues 
of Mcst and the elements Uij of the orthonormal matrix U are given by 
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with i, j = 1, . . . ,n. In this particular example we have 



fij = uj + 2c — 2c cos 



n+1 



uj'^ + 4csin^ 



jvr 



2(n + i; 



as eigenvalues of Acst- The operators are constructed as shown above and the entire Wigner 
quantization procedure can be reproduced with the given information. 

Kraw^tchouk interaction [2] . Another example of a matrix for which the eigenvalues and 

eigenvectors have analytically closed expressions, is the Krawtchouk matrix, given by 



Mk 



-El 




-El 
Fi 

—E2 







— E2 

F2 



-E, 



-En-l 
Fn-1 J 



where 



Fr = -p)yjr{n - r), = {n - l)p + (1 - 2p)r, 

for r = 0, . . . , n — 1. The parameter p lies between and 1. The spectral decomposition of such 
a matrix is known. The eigenvalues can be written as = j — 1, for j = 1, . . . , n, and, following 
the notation of Section [2l the element at position = l,...,re) of the orthonormal 

matrix U is given by Ki^i[j — 1), with 



K^{3) 



n — 1 
i 



n — 1 
j 



P 



n—i—j — l 



1/2 min(ij) /jx /jx 



k=0 



These elements Ki{j) are evaluations of normalized Krawtchouk polynomials, which explains 
the name of the matrix. 

Consider the system with 



A = io^I + cMk. 

Then this Hamiltonian is analytically solvable and can be rewritten as 



(13) 



r=l 



vl mu)'^ + cmpln — 1) 9 cm(l — 2p) , 

+ -ll + ^ rql 

2m 2 2 



cra\Jp(\ — p) 



^Jr{n - r){qrqr+i + qr+iQr)- 



r=l 



It is clear that this simplifies a lot if the parameter p is chosen to be ^ . The Hamiltonian then 
takes the form that is presented in the introductory section. 

Following the general procedure of Section [2l we can rewrite the problem as (jlip with 



fij = (j^ + cXj = uj"^ + c{j — 1), 



1,2, 



,n. 



The Krawtchouk interaction matrix (113p is clearly real and symmetric. Moreover, the quantities 
are well defined for all j = 1, . . . , n and the matrix (jlSp is positive definite whenever uj ^ 0. 
Other types of interaction matrices connected to discrete orthogonal polynomials can be 
found in [2] . We note that the interaction matrices in the previous examples can be written as 



A = io^I + cM. 



(14) 
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From now on we will always work with an interaction matrix that can be written in the form (jl4p . 
As a consequence, whenever we adopt results like Q or (fTUj) from Section [51 one must first set 
fjLj = u;^ + cXj. The entire Wigner quantization procedure described in this section can be 
reproduced using this transformation. 

It is now the task to find operator solutions for that satisfy the equation 



= ±2v'w2 + cAfe a^, k = l,2,...,n. (15) 

As stated above, we are able to express the wanted operators in terms of generators of Lie super- 
algebras. We can then find solutions in specific representation spaces of these Lie superalgebras. 



E 



+ cXj{aj,aj },< 



4 Lie superalgebra solutions 
4.1 The £(l(l|n) solution 

The Lie superalgebra 5[(l|n) has basis elements ejk with j,k = 0,1, ... , n. The solutions we are 
about to find, will be in terms of the odd elements ejo and cqj with j = 1, . . . ,n. The remaining 
basis elements are called even elements. The odd elements have degree 1, the even elements 
have degree 0. The commutation and anti-commutation relations in g[(l|n) are determined by 
the Lie superalgebra bracket 

le^„ekl} = 5,keu - (-l)'i^s(^-)'^^g(^'=')5,;efc,-. 

In a representation of this Lie superalgebra, the bracket [[x,y] stands for an anti-commutator 
if X and y are both odd elements of 5((l|n), and for a commutator otherwise. We can use a star 
condition for gl(l|n) that is fixed by a signature a = (ui, . . . ,cr.«), a sequence of plus or minus 
signs, and by 

{eoj)^ = ajCjo, j = l,...,n. (16) 

We will restrict ourselves to the case where all cij's are equal to -|-1, since this corresponds to 
the real form u(l|n). In this case it is known that finite-dimensional unitary representations 
exist [T5] . 

Solutions of (jl5|) in terms of generators of 3[(l|n) are known. They have been constructed 
for a fixed interaction matrix in [9]. Therefore, it is not necessary to copy the entire analysis 
here. The solutions are of the form 



a, = J e,o, a, = sign(/?,) W eo„ J = l,...,n, (17) 



where the /3j's are given by (n > 1) 

^ n 

fij = -^Cj2 + cXj + ^—^ + ""^k- (18) 

k=l 

It is straightforward to verify that the Hamiltonian ([9|) can be written as 

H = hheoo + J2PjejA , (19) 
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where we have used the notation (3 = Pk to indicate that this is a constant. We want 

aU fjj's in equation (|16|) to be equal to 1, which is, together with the adjointness condition 



aJ, equivalent to saying that the values (3j need to be positive. Examining the form 



of f3j we see that it is equal to — ywM-cAj plus some average value of the roots of the quantities 
(u;^ + cXj). Thus, it seems reasonable that half of the /3j's will be positive and half of them will 
be negative. However, it is possible to prove that this is not always the case. More concretely, 
we will need to assume that the coupling strength c is small enough. We shall refer to this as 
"weak coupling". 

Note that the condition that all Pj's are positive does not stem from the algebraic solution 
of (lisp by means of (I17p . but only from the requirement of the star condition (eoj)^ = ejo, since 
we are primarily interested in unitary representations of the real form u(l|n) of 0l(l|n). 

Krawtchouk interaction. Remember that we are considering a system of n harmonic 
oscillators that are coupled by a certain interaction matrix. Finding operator solutions for the 
Hamiltonian of this system treated as a WQS was proved to be equivalent to finding operators 
that satisfy the relations (jl5p . This equation contains the eigenvalues of A and the operators 
are dependent on the eigenvectors of A. In the specific case of Krawtchouk interaction, we know 
that the jth eigenvalue Aj is equal to j — 1, with j = 1,2, ... ,n. We will find an upper bound 
for the coupling strength c so that in this case all the /3j's are positive. 

In order to find an upper bound for the value of c, we need the following property. 

Lemma 1. For C > ^"-^g^ , we have the inequality 



Y.VcT]>{n + l)Jc + --l. 
j=o ^ 

Note that C denotes an arbitrary positive constant in this lemma, it is not the coupling strength. 

Proof. The proof goes by induction on n. For n = 1 and n = 2 the property is trivial, as one 
sees for example from 

Vc + VcTi + VcT2 > sVc. 

For larger n, we first notice that 



Vc + VcT^>2jc+'^-i 

if C > — 1)^. This can be verified by solving the previous inequality to C as if it were an 
equality. Taking the square of both sides twice results in the given boundary for C. Consequently, 
one finds 

n n—2 

0=0 i=o 



>2^/C + ^-l + (n-l)^(C+l) + ^-l = (n + l)yC7+^-l, 

where induction is used to justify the inequality. This is possible because if C> then 
surely C + 1 > ^"^g ^ as long as n > 1. ■ 

It is then possible to construct an upper bound for the coupling strength c, as is shown in 
the following proposition. 
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Proposition 1. Assume that the eigenvalues of the interaction matrix A are equal to fij = 
uj^ + cXj, with Xj = j — 1 (j = 1, 2, . . . , n). An upper bound for the coupling strength c is then 
given by 

2(2n - 3)^2 
^ (n - l)(n2 - 3n + 4)' 

If c satisfies this condition, then all the f3j 's given in equation (jlSp are positive. 

Proof. First of all, since l3j — jSj-i = ^/lo"^ + c{j — 2) — y^w^ + c (j — 1) < 0, the row j3j 
(j = 1, . . . , n) is decreasing. All of the /3j's will thus be positive if and only if /3„ is positive. An 
equivalent condition is determined by 

n~l I 2 I 2 

Pn>0 ^ I^Y^ + i > (^-1)Y^ + ^-1- (20) 

j=0 

To prove this inequality, we want to use Lemma [1] for n — 2. Therefore, we need to check if the 
condition of the lemma is satisfied. For ?i > 2 we have that 

^2 (n- l)(n2 - 3?i + 4) (n-6)2 
~ ^ 2(2n - 3) - 16 ■ 

Since we will only consider systems with at least two coupled harmonic oscillators (n > 2), 
Lemma [1] is applicable: 

^-+,>(n-l)^- + --2 + ^- + n-L 

j=0 

By demanding that the r.h.s. of the previous equation is larger than or equal to {n—\)\J ^+n— 1, 
we ensure that f3n is positive. A simple calculation shows that this is true for values of c that 
are smaller than or equal to the upper bound given in this proposition. ■ 

The upper bound for c/uJ^ is of the order A/n^. An idea of how accurate our approximation of 
the boundary value is, can be found in Table [TJ In this table, c„ denotes the highest value for the 
coupling strength c for which /3„ and hence all the /3j's are positive. These values can be found 
by solving equation (j20p numerically to c. The boundary value as proposed in Proposition [T] is 
denoted by c^. 



Table 1. Critical values CnjuP' in the case \j = j — 1 



n 




Cn 


Cn 
Cn 


n 


Cn 


Cn 


Cn 
Cn 


4 


0.41667 


1.27357 


0.32717 


9 


0.06466 


0.08639 


0.74843 


5 


0.25000 


0.51723 


0.48334 


10 


0.05105 


0.06562 


0.77802 


6 


0.16364 


0.27857 


0.58742 


20 


0.01132 


0.01259 


0.89893 


7 


0.11458 


0.17391 


0.65886 


50 


0.00168 


0.00175 


0.96186 


8 


0.08442 


0.11887 


0.71013 


100 


0.00041 


0.00042 


0.98130 



For example, if n = 8, all /3j's will be positive if the coupling strength c < 0. 11887 w^. The 
boundary value from Proposition [T] is a little more pessimistic (c < 0.08442 cj^), but not too 
much. One might believe from these observations that the fraction c„/c„ tends to 1 for large n. 
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Constant interaction. The problem of finding a boundary value for the coupling strength c 
so that all the /?j's are positive was discussed for constant interaction in [TO]. The authors propose 
an estimation of the upper bound for the coupling strength. Moreover, for n = 4, . . . , 21 they 
give exact values of this upper bound [TOl Table 1, page 22]. 

4.2 The osi3(l|2n) solution 

Apart from the solution in terms of generators of the Lie superalgebra 3((l|n), we can also 
express a class of solutions of (jl5p by means of 05p(l|2?i) generators. It is known (see [16]) that 
this Lie superalgebra is generated by a set of 2n paraboson operators {j = 1,2, ... ,n) that 
satisfy the relations 

[{bl, bl}, bt] ={e- OS.ibl + (6 - v)6,4 (21) 

In these triple relations, j, k and / are elements from the set {1, 2, . . . ,n} and rj,^,e £ {+, — } (to 
be interpreted as +1 and —1 in the algebraic expressions (e — C) and {e — t])). The elements b^ 
are the odd elements of the superalgebra, while the even elements are formed by taking anti- 
commutators {b^j,b^}. 

Using the relations (j21|) it is then easy to check (see also [11]) that the operators 

with j = 1,2, ... ,n indeed satisfy equation (fT5]) . The Hamiltonian ([9]) then takes the following 
form: 

n , n 

H = ^-,f^7^,{al,aj} = hY,\f^^+^jhj, (22) 
j=i j=i 

where we have introduced the notation hj = {a^,a~}/2 = {6^,6j}/2. The Cartan subalgebra 
of osp(l|2n) is spanned by the n elements hj {j = 1,2, . . . ,n). 

5 The spectrum of in a class of representations 

In order to study the spectrum of the Hamiltonian H in terms of the g[(l|n) or osp(l|2n) 
solutions, it is necessary to work with specific representations of those Lie superalgebras. An 
approach to this problem with respect to the Fock-type representations W{p) of 0l(l|?i) was 
given in [9]. Here, we will work with another type of representations V{p). 

5.1 The gl(l|n) representations V{p) 

Before analyzing the spectrum of H, we will summarize the main features of the representa- 
tions V{p) of g[(l|n). First of all, they are finite-dimensional, unitary representations. For any 
natural number p, the basis vectors of V{p) are given by |18j 

v{e;r) =v{e;ri,r2,. . . ,rn), 

with 6 G {0, 1}, Ti G {0, 1,2, . . .} and 6* + ri + • • • + r„ = p- The dimension of the vector 
space V{p) equals 

fp + n—l\ /p + n — 2 
\ n — 1 J \ n — 1 
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in which both terms represent the number of basis vectors for 6 = and 6 = 1 respectively. The 
action of the 0l(l|n) generators on these basis vectors can be determined [T5j. We will give the 
actions of the diagonal elements eoo and e^fc since only these actions will be needed to find the 
spectrum of the Hamiltonian: 

eoov{6;r) = 6v{6;r), 
ekkv{6;r) = rkv{6;r). 

In terms of the 5[(l|n) generators, the Hamiltonian takes the form (jl9p : 

H = h ^/3eoo + f^j^h 

Clearly, looking at the actions of the elements eoo and e^fc, the vectors v{6;r) are eigenvectors 
for H: 

Hv{6;r) = hE^v{6]r), 

where the eigenvalues hEj. are determined by 

n 



Er = pp-2_^ ^Lo^ + cXj rj . (23) 
i=i 

This can be established by noting that 

n -. n 

and by using the fact that + ri + • • • + r„ = p. 

In the case where there is no coupling (c = 0), all the /3j's become the same. It follows that 
Pj = uj/{n — 1) and P = nuj/{n — 1). In this case, we thus see that the eigenvalues of H are 

\n-l 

So in fact, there are two eigenvalues. The lowest one, for 6* = 0, has multiphcity {^\'!:\^). The 
highest eigenvalue has multiplicity (^^"Y^) • 

Our main interest lies in the weak coupling case, where < c < c„. The energy levels are 
easily computed through equation ([23]) . The result for n = 4,p = 2 and uj = h = 1 can be 
seen in Fig. [H where we have chosen to compare the systems with constant and Krawtchouk 
interaction. 

Both figures look quite similar, but there are some differences. We see that in general all 
eigenvalues are different, but for specific values of c some of the energy levels cross each other. 
For these values of c, the multiplicity of some of the eigenvalues is higher than 1. In the constant 
interaction case, we see that energy levels can cross if we restrict ourselves to, say, = 0. Also, 
note that there are indeed only two eigenvalues in the case without coupling (c = 0). 

The figure also suggests that the lowest energy level tends to zero as the coupling strength 
reaches Cn- In order to prove this, we need to know the lowest energy level. First, we note 
that j3n < (3j for all j as soon as A„ > \j for all j. We can always choose A„ to be the largest 
eigenvalue, so we can assume that /3„ is the smallest of all (ij. Next, the formula 

n 

Er = fie + Y,Pjrj 
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Constant interaction Krawtchouk interaction 




Figure 1. Spectrum of the Hamiltonian (fT9|) in the 0[(l|n) representation V{p) for n = A, p = 2 and 
Lu = h = 1, as a function of the couphng constant c. The left figure belongs to the system with constant 
interaction, where the Xj are given by equation (|12p . The right figure represents the Krawtchouk case, 
with Xj — j — 1. 



gives the energy levels as a sum of p terms. Thus, the lowest energy level arises when all p 
terms are equal to The definition of Cn tells us that /3„ = for c = c^, so the lowest energy 
level pPn tends to zero when c approaches c^i. 



5.2 The 03^3 (1 1 2n) representations V{p) 

We will also take a look at the (infinite-dimensional) representations V{p) of 05p(l|2n), with 
lowest weight (|, . . . , |). Such a representation is a unitary, irreducible representation (unirrep) 
if and only if p G {1, 2, . . . , n — 1} or p > n — 1 \n\ Theorem 7]. In literature, where 05p(l|2?i) 
is related to the n-paraboson algebra, the parameter p is sometimes referred to as the order of 
the parastatistics. A basis for the representations V{p) was given in [T7], and consists of all 
Gelfand-Zetlin patterns for partitions of length at most n. These GZ-vectors have the following 
form: 



\m) = |m)" 



min 

"il,n-l 

mil 



\m\ 



\m 



in-1 



The top line of the GZ-vectors is any partition into at most p parts, where p is the label of the 
representation. This partition is denoted by [m]". For basic information about partitions, we 
refer to [19]. The other elements of the GZ-vectors, denoted by \m)^~^ satisfy the betweenness 
conditions 



l<i<j<n — 1. 



The actions of the 05p(l|2?i) generators on these basis vectors are known [17]. In particular, the 
action of the diagonal elements hj is given by 



hj\m) 



P 



+ rurj — m^j-i j |m). 



The Hamiltonian of the system in terms of 05p(l|2n) generators was given by equation ([22]): 



H = h Y + cXj hj 
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From the action of the diagonal elements hj it is clear that the vectors \m) are eigenvectors of 
the Hamiltonian. We can write 

H\m) = hEm\m), 

in which Em stands for 



Em = '^\J^'^ + C\j ( 2 + ~ "T-rj-l I • (24) 

j=l \ r=l r=l / 

In the case without coupling (c = 0) we see that the eigenvalues simplify significantly and they 
can be written in the form 



n 
r=l 



The summation in this expression is in fact the weight of the partition [m]". This weight can 

(p) 

be any positive integer k, which we shall call the height of the eigenvalue Ef^ . This means that 
for c = there is an infinite amount of eigenvalues, that can be written as 



E^P^ = hu;(!^ + k), /fc = 0,l,2,.... 

The multiplicity fj.{E^^) of each eigenvalue can be determined with the help of some theoretical 

arguments. First of all, fi{E^^) will be equal to the total number of GZ-vectors with a partition u 
in the top row, where v is any partition of k into at most p parts. Let v' be the conjugate partition 
of v [19]. It is known (see for example [201 Section 4.6]) that the representation of gl{n) that 
is labelled by the partition u has dimension (^,) , where we have used the generalization of the 
binomial coefficient for a partition |19[ page 45]. This is defined by 

_ -Q X-C{i,j) 



where c{i,j) = j —i and h{i,j) = Vi + u'^ —i — j are the content and the hook length of (i, j) 
respectively. So for a given partition z^, the number of GZ-patterns that have u in the top row 
equals (^,) . This implies that the multiplicity of each eigenvalue is equal to 



u, \u\=k,l{u)<\p\ 



The ceiling function \p\ is used to cover the cases where n — \ < p < n. So we have found 
that the energy levels for c = are equidistant with spacing huj and the multiplicities of the 
eigenvalues can be computed through equation ([25]) . 

In the case with actual coupling (c ^ 0) the eigenvalues can be found by equation (|24p . 
Unlike the weak coupling case in the representations V{p) of gl(l|n), the multiplicities of the 
eigenvalues are not all equal to one. Any two basis vectors |m) and |m') that are subject to 

3 j 

'^nirj = ^m',.j Vj = l,2,...,n (26) 

r=l r=l 
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yield the same eigenvalue hEm- So the multiplicity of an eigenvalue hEm is equal to the number 
of basis vectors for which the sum of the elements on row j is equal to the sum of the elements 
on row j of |m), for every j. For example, the vectors 



\m) 



5 \ 
4 
2 



1 



and 



\m = 



3 2 \ 
3 1 
2 



1 



yield the same eigenvalue. From this it is also clear that the total number of distinct eigenvalues 
at height k is equal to that number for p = 1. Indeed, every vector with p > 1 can be associated 
with a vector with p = 1 for which equation (|26p holds, as can be seen in the previous example. 
Moreover, all eigenvalues in the case p = 1 clearly have multiplicity 1 (in the generic case when 
all \j are distinct). We know what the number of eigenvalues at height k for p = 1 is, namely 



E 

I/, \i'\=k,l{i>)<l 



n 



n 

{ky 



k 

n 



n + j — 1 
j 



where the default value for A; = is equal to 1. The latter product is nothing more than the 
binomial coefficient ("^^^"^), which shows that it is an integer. 

It is now clear that some eigenvalues have multiplicity greater than 1. Furthermore it is 
possible that some of the energy levels cross each other, just as in the g[(l|n) case. This means 
that for specific values of c there are some eigenvalues for which the multiplicity is even higher. 
It would be inappropriate to try to compute these values of c. Let us instead look at Fig. [21 
where we have plotted a part of the energy spectrum for n = 4,p = 2, u; = ^= l and \j = j — 1, 
to visualize things. Recall that we are dealing with an infinite spectrum. Therefore, we will 
only plot the spectrum up to height k, for k = \ and k = 2. 



k = 1 



Figure 2. Spectrum of the Hamiltonian ((22)) in the osp(l|2n) representation V{p) ior n — A, p ~ 2, 
Lu = h = I and Xj = j — 1, as a function of the coupling constant c. The left figure gives the spectrum 
up to height fc = 1, the image on the right goes one step higher (fc = 2). The total spectrum is infinite. 

The eigenvalues on height and 1 all have multiplicity 1 for c > and they never cross. 
The figure on the right shows the energy values for = 2 as well, where we both have higher 
multiplicities and crossing energy levels. Six of the ten distinct energy levels at height 2 have 
multiplicity 2. 

6 Relation to canonical quantization 

Having treated the quantum mechanical system ([1]) as a Wigner quantum system in general, 
it is interesting to consider the canonical case as a special case. In this case, one assumes the 
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canonical commutation relations (CCRs) 

[qr,qs]=0, [Pr,Ps]=0, [qr,Ps] = ih5rs- (27) 

Similarly as before, we write A = U^DU, where U is the orthonormal matrix with the eigenvec- 
tors of A as rows and D is the diagonal matrix with the eigenvalues of A on the main diagonal. 
We define new operators Q = U{qi . . . Qn)'^ and P = U{pi . . .pn)^ ■ These are subject to the 
same commutation relations as in equations (j27p and yield a new expression for the Hamiltonian: 



We then take the linear combinations ^ of the operators Q and -P, and thus create new opera- 
tors . The Hamiltonian can be written in terms of these new operators: 

" h 



i=i 

Using the canonical commutation relations of Q and P, one finds that the operators satisfy 
the usual boson commutation relations: 

As before, it can also be verified that 



[H, af] = ±h^LO^ + cXj af, j = 1, 2, . . . , n. (28) 

These are in fact the compatibility conditions derived in the general case, interpreting the system 
as a WQS. These CCs are also valid in the canonical case, since the CCRs imply the CCs. 

We will now define the n-boson Fock space, which is equivalent to the representation V{1) 
of osp(l|2?i). Since p = 1 represents the canonical case, we will find a correspondence between 
the basis vectors of V{1) and the basis vectors of the n-boson Fock space. The latter are 
constructed from a vacuum vector |0), with 

(0|0) = 1, aj\0) = 0. 

The other (orthogonal and normalized) basis vectors are then defined by 

\ku...,K)=^^^^^==^$p\0), (29) 

where ki,...,kn G Z+. We need to find a correspondence between these vectors and the 
Gelfand-Zetlin basis vectors of the representation V{1) of osp(l|2n), generally denoted by 

m„ • • • \ 
m„_i • • • 



1712 

mi 



(30) 



Using [a^. ,a^] = 1, one finds that 



n 



2 
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where we have used the notation Eq to indicate the lowest energy level. This is the lowest 
energy state of the system. The higher energy levels can be calculated using equation (|28|) in 
a straightforward way. This results in 



H\ki, ...,kn) = h^EQ + ^/cj^w^ + cXj j \ki, . . . , A;.„). 



By comparison with equation (j24p for p = I, one finds that 

j i-i 
kj = TTirj — rrirj-i = rrij — rrij^i. 

r=l r=l 

Thus we have: 

Proposition 2. The n-boson Fock space and the osp(l|277,) representation space V{1) are equiva- 
lent and their basis vectors ()29p and ()30p are related by kj = mj — mj-i. 



7 Conclusion 

To conclude, we have in this paper considered the quantization of a system of harmonic oscillators 
as a Wigner quantum system. The quadratic coupling terms have been characterized by an 
interaction matrix. For such systems, the Wigner quantization procedure can be performed 
completely (Theorem [T|), leading to a set of algebraic triple relations (jlip as compatibility 
conditions. These relations have particular solutions in terms of generators of the Lie superalgeb- 
ras 0[(l|n) or osp(l|2n). Then the unitary representations of these Lie superalgebras play an 
important role: the algebraic generators, and thus also the physical operators corresponding to 
observables, act in these representations. For some classes of representations, the spectrum of 
the Hamiltonian operator is determined explicitly, and discussed. 

As leading examples throughout the paper, we consider two analytically solvable systems. The 
first is a classical one, describing a linear chain of harmonic oscillators coupled by a harmonic 
nearest-neighbour interaction. The second is a relatively new system, again describing a chain of 
harmonic oscillators, but this time the nearest-neighbour coupling is a "Krawtchouk coupling" . 

The original results in this paper are: the proof that all Hamiltonians with quadratic interac- 
tion terms can be reduced to an algebraic set of relations under Wigner quantization (Section [2]); 
the conditions for the coupling strength in the case of Krawtchouk interaction for the g[(l|n) 
solution (Section |4]); the determination and discussion of the energy spectrum for specific rep- 
resentations of the Lie superalgebras 0l(l|n) and osp(l|2n) (Section [5]). Wigner quantization is 
an extension of canonical quantization, so the canonical case appears as one particular solution 
of the various solutions allowed by Wigner quantization. For the systems under consideration 
here, the canonical case corresponds to the representation V{1) of the 05p(l|2n) solution. This 
correspondence is explained in Section [6l 

In this paper we have dealt with the Wigner quantization of Hamiltonians of the form ([3]), 
where the interaction terms are explicitly present, so there is 'dynamical interaction'. It is 
interesting to note that Wigner quantization also allows another peculiar feature. In [21], one 
considers the Wigner quantization of a system of free oscillators (i.e. the Hamiltonian contains 
no dynamical interaction terms), and it is shown that certain solutions of the compatibility 
conditions can be interpreted as being responsible for statistical interactions in the system. So 
in such a case the purely algebraic quantization condition implies a statistical interaction in 
a free system. 
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